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Abstract 

The developments toward a general nultilevel 
optimization capability and results for a three- 
level structural optimization are described. The 
latter is considered a major stage in the method 
development because the addition of more levels 
beyond three does not Introduce any new qualitat- 
ive elements, so that a three-level implemen- 
tation is qualitatively equivalent to a multi- 
level implementation. 

The method partitions a structure into a 
number of substructuring levels where each sub- 
structure corresponds to a subsystem in the 
general case of an engineering system. The 
method is Illustrated by a portal framework that 
decomposes into individual beams. Each beam is a 
box that can be further decomposed into stiffened 
plates. Consequently, substructuring for this 
example spans three different levels: the bottom 

level of finite elements representing the plates, 
an intermediate level of beams treated as sub- 
structures, and the top level, for the assembled 
structure. This example is an extension of a 
case presented previously which was limited to 
two levels. Further extensions would add only 
more intermediate substructuring levels; there- 
fore, the three-level case is qualitatively 
complete. 

Nomenclature 

Quantities 

A Cross-sectional area 

C Cumulative constraint (equation 11) 

c Capacity: limitation on the ability 

to meet a particular demand D (e.g., 
allowable stress) 

D Demand: a physical quantity the 

structure is required to have, to 
support, or to be subjected to in 
order to perform its function (e.g., 
stress) 
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F Objective function 

f( ) Functional relation 

g Vector of constraint functions, gj 

h Equality constraints defined by 

equation 10 

I Cross-sectional moment of inertia 

K Stiffness matrix 

KS A function defined by equation 11 

L Lower bound on X t including move 

1 imits 

M Mass (or "middle" when superscript) 

Q Boundary forces on SS 

q Number of the diagonal and off- 

diagonal symmetric entries in K 

SS A substructure (including the extremes 

of the assembled structure and a 

single structural element) 

STOC Acronym: subject to constraints 

U Upper bound on X including move limits 

X Vector of design variables, X* 

Y Vector* containing those entries of K 

and the mass M that are to be held 

constant in an SS optimization 

n A vector defined by equation 15 

p A user-controlled constant in the KS 

function 

a Increment of a variable (see 

definition of subscript o) 

Indices, Subscripts, and Superscripts 

Overbar Denotes an optimal quantity 

B Superscript for bottom level 

b Superscript to denote an association 

with the SS boundary 

e Subscript to identify an extrapolated 

value 
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1 Index or superscript to identify the 

substructure level, i = 1 to i max , 
see Figure 1 

j Index or superscript to identify the 

substructure position, counting from 
left, at level i, see Figure 1, or 
subscript of an entry in g in optimi- 
zation without decomposition 

k Index or superscript to identify index 

j of the parent SS, see Figure 1 

1 Equivalent of k, see Figure 1 

M Superscript for middle level 

m Equivalent of i , m = i - 1 , see Figure 1 

n Equivalent of i, n = m-1, see Figure 1 

o Subscript to identify an- original 

value from which a is measured 

p Equivalent of k, see Figure 1 

r Subscript of an entry in Q, r = 1 to R 

s Subscript of an entry in h, s = 1 to S 

T As superscript denotes a transposed 

vector 

t Subscript of an entry in X, t * 1 to T 

w Subscript of an entry in g in multi- 

level optimization, w = 1 to W 

z Subscript of an entry in Y, z = 1 to Z 

Introduction 

Formal optimization methods applied to 
realistic structures built-up of many components 
and carrying a large number of loading cases are 
hindered by an excessive number of design vari- 
ables and constraints. They frequently become 
too costly and unmanageable, and can easily 
saturate even the largest computers available 
today or in the foreseeable future. An obvious 
remedy is to break the large optimization problem 
into several smaller subproblems and a coordi- 
nation problem formulated to preserve couplings 
among these subproblems. A very Important ben- 
efit of such an approach, in addition to making 
the whole problem more tractable, is preservation 
of the customary organization of the design 
office in which many engineers work concurrently 
on different parts of the problem. Therefore, 

research has been directed toward multilevel 
optimization methods that decompose large 
problems into a hierarchically related set of 
smaller subproblems while preserving their 
coupling. This approach meshes well with the 
recent trends in computer technology toward 
computing distributed over a network of computers 
whose characteristics may be matched to individ- 
ual subproblems for more efficiency and con- 
venience. Moreover, the decomposition approach 
is natural in an engineering organization. Since 
engineers tend to cluster into groups concentrat- 
ing on parts of a project in order to develop 
broad work front to shorten the development time. 


A number of procedures for implementation of 
the foregoing approach has been proposed for 
structural applications, (e.g., refs. 1, 2, and 
3). A multilevel optimization with decomposition 
has also been formulated in a general manner for 
use in engineering system design (ref. 4) con- 
cerned with the "trade-offs" among various physi- 
cal subsystems that may be governed Dy different 
engineering disciplines. The* unique feature of 
the algorithm proposed in reference 4 is the use 
of the optimum sensitivity derivatives introduced 
in reference 5 as means to approximate the coup- 
lings among the subsystems. 

When the system optimization formulation 
established in reference 4 is applied to struc- 
tural optimization, it's analysis part coincides 
with a general, multilevel substructuring (e.g., 
refs. 6, 7, and 8). In the simplest case, the 
system acquires a meaning of a complete structure 
and each subsystem corresponds to a single struc- 
tural component that may be represented by a 
single finite element. This is a two-level, 
structural optimization whose algorithm was 
illustrated by an example of a framework reported 
in reference 9. This served as a verification of 
the general purpose algorithm laid out in refer- 
ence 4. 

Since the general algorithm presented in 
reference 4 allows a theoretically unlimited 
number of hierarchical subsystem levels in the 
decomposition, its continuing development 
requires verification by applications of more 
than two levels. The purpose of this paper is to 
report such a verification using structural opti- 
mization with a three-level decomposition. In 
this decomposition, the highest level corresponds 
to the assembled structure, the level below 
corresponds to substructures, and the third, 
bottom level, represents the structural 
components that make up the substructures. 

Extension of the scheme beyond the three 
levels would require more levels of nested sub- 
structures sandwiched between the top and bottom 
levels. Thus, such an extension would not add 
any qualitatively new subsystems to the scheme 
and, therefore, one may regard the three-level 
optimization as the simplest case of the most 
general multilevel optimization. 

The paper presents an multilevel algorithm 
for structural optimization and its verification 
by a three-level optimization of a framework 
structure. 

One-Level Optimization 

An optimization formulation without decompo- 
sition serves as a reference from which the 
multilevel optimization algorithm is derived. 

The optimization is defined in terms of; 
the design variables, X^, which are the cross- 
sectional dimensions of the structural compo- 
nents; the objective function F { X) that can be 
any computable function of these variables 
(structural mass is the frequent choice); and the 
constraints, gj(X), Imposed on the the behavior 
variables to account for the potential failure 
modes. It is useful to distinguish here between 
the "local" constraints such as local buckling 
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that depend predominantly on the component behav- 
ior and the "global" constraints that primarily 
depend on the characteristics of the assembled 
structure (e.g., displacements or overall buck- 
ling). However, in the one-level optimization 
formulation both constraint categories will be 
treated in the same way. 

Writing the constraint functions as 

g = (D/c - 1) _< 0 (jj 

the optimization proolem in a standard formu- 
lation is 


where K bi+ ^>J denotes the boundary stiffness 
matrix of the jth daughter substructure. 

The stiffness matrix . of .the SSijk is 
assembled of the matrices K bl+1 »J - the fact 
symbolically expressed as 

K ij = S (K bi+ l.j) - (4) 

where S stands for an appropriate stiffness 
summation operator. 

Similarly, the mass, M^, of SSijk is a 
simple sum of the daughter masses 


find min F(X) 
X 


M 1 -} = l 


STOC gj(X) £ 0 (2) 

and requires a search of the n-dimensional design 
space considering all the design variables and 
constraints concurrently. In contrast, an algor- 
ithm presented in the next section breaks the 
problem into a number of search operations, each 
concerned with a smaller number of design vari- 
ables and constraints. 


For a substructure SSijk that is at the ultimate 
level of detail and, therefore, is pot further 
subdivide, the stiffness matrix, K b1 J, and the 
mass, M n J, derive directly frqm the cross- 
sectional dimensions denoted by X^J so that 


K bij 8 K ij . f(xij) 


( 6 ) 


A Multilevel Optimization 
Preliminary Definitions 


and 


Mij = f(xiJ) 


( 7 ) 


The diagram in Figure 1 shows a structure 
decomposed into several levels of substructures. 
The term "substructure" will refer to any entity 
in this decomposition scheme other than the full, 
assembled structure represented by the box on the 
top of the pyramid. In the limit, then, a sub- 
structure may be a single structural component 
representing the ultimate geometrical detail 
appropriate to the problem at hand. The sub- 
structure levels are numbered from 1 on the top 
to imax at the bottom. The hierarchical nature 
of the scheme instigates the use of a term 
"parent" to the structure at level i which, in 
turn, is decomposed into a number of "daughter" 
substructures at level 1+1. A daughter may have 
only one parent and that parent must be at the 
level immediately above. Thus, it will be con- 
venient to label each substructure SSijk, where i 
denotes the level, j defines the position at the 
level i counting from the left, and k identifies 
the parent's position at the level 1-1. The sub- 
structure occupying the last position in a 
particular parent-daughter succession represent 
the ultimate level of detail at which the 
decomposition stops. There Is no requirement 

that all such structures must be at the same 
bottom level imax. Any particular parent- 

daughter succession line may end at 1 < i <= 
imax. 

Substructuring analysis (e.g., refs. 6, 7, 
8) of a parent’' substructure SSijk, yields the 
forces acting on the boundary of each daughter 
substructure. In the general case these forces, 
collected in a vector Q, depend on the stiffness 
properties of all the daughter substructures so 
that 

Q = f(K bi+ l.J) 


The same cross-sectional dimensions and other 
appropriate geometrical and material data can be 
entered together with the forces Q into a 
procedure to calculate stresses, internal forces, 
and the critical stresses (and/or internal 
forces) for local buckling. This information, 
together with the displacements calculated in the 
analysis of the entire structure and each sub- 
structure and overall elastic stability analysis 
results, describe the static elastic behavior 
that may be subjected to constraints in the opti- 
mi zation. 

Although the foregoing definition of sub- 
structuring analysis is based on the stiffness 
approach, the use of a finite element analysis is 
not mandatory for the multilevel optimization 
algorithm whose description will follow. As far 
as that algorithm is concerned, the analysis is a 
"black box" where only the inputs and outputs are 
important but not the content. 

Multilevel Optimization Algorithm 

With the substructuring scheme and analysis 
established in the foregoing, this section 
describes the optimization algorithm itself. The 
essentials of the computer implementation are 
also given. 

Basic Concept 

The basic concept of the multilevel struc- 
tural optimization introduced in reference 9 is 
based on the well-known property of a substruc- 
ture that the elements of the boundary stiffness 
matrices of its daughters on the right-hand side 
in equation 4 can be changed in such a way that 
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its stiffness matrix, K b ’J, will be held con- 
stant and, therefore, the boundary forces, Q, 
will remain unchanged. In the case of a sub- 
structure that is not a parent, the same applies 
to its cross-sectional dimensions and its stiff- 
ness matrix, K’j, provided, of course, that the 
number of symmetric entries, q’j, in K 1 J is 
smaller than the number of cross-sectional dimen- 
sions, T’J, in the vector X 1 -!. 

qiO < lU ,o, 


When this inequality holds, there is a design 
freedom to "tailor" the substructure to a set of 
prescribed stiffness properties. Otherwise, if 


qij > T’j 


(9) 


0^ 


j’J 


hiJ 


Vector of the boundary forces, 
Qij; r = 1 through R’j. 


•Vector of the inequality con- 
straint functions, g 1 ^, w = 1 

through W 1J . 

Vector of the equality con- 
straint functions, h 1 J ; s = 

1 through si j . . These . con- 
straints link Y 1 j and X’j so 
that considering equation 6 


hi j 
hij 


yij . f(icbU) = o a) 

Y 1 J - f ( X i j ) = 0 b) 


then the set of prescribed stiffness properties 
either defines the cross-sectional dimensions 
uniquely or cannot be physically realized. 

Based on the property discussed in the fore- 
going, the optimization scheme introduced in 
reference 9 and extended here uses the. elements 
of the boundary stiffness matrices, K b ’ + l»j, of 
the daughters as design variables of their parent 
substructure. That definition of the design 
variables applies to all substructures except 
those that have no daughters. Their design vari- 
ables are their cross-sectional dimensions. 

Optimization At The Most Detailed Level 

Introduction of the optimization algorithm 
begins at the level of the most detailed sub- 
structures that are not further subdivided. It 
is assumed that a complete, top-down, substruc- 
turing analysis has been carried out so that for 
an. SSijk one has computed its boundary forces 
Q’J. The prerequisite completion of the 
analysis implies that the entire structure has 
been initialized. Consequently, the SSijk has 
its mass, M' J , cross-sectional dimensions, 
X’J, and the entries of its stiffness matrix, 

K 1 j , as the given quantities. 

The optimization problem to be solved for 
SSijk calls for the cross-sectional dimensions to 
be treated as design variables and manipulated so 
as to minimize a measure of the constraint viol- 
ations local .t.o SSijk while holding the elements 
of matrix K 1 J and the mass M’J constant and 
staying within the side constraints. 

To formalize this, define: 


L^, 


where f( ) denotes a known, 
computable function. 


U 1 j 


Vectors of lower 
limits on the design 

x’J. 


and upper 
variables 


C 1J Cumulative constraint, .a. single 

valued function of g 1 j , con- 
tinuous and differentiable, 
having the property of being 
positive wjipn at least one con- 
straint gJJ is positive (that 

is violated in the convention 
adopted here). Following refer- 
ence 9 the cumulative constraint 
is chosen in form of the 
Kresselmeir-Steinhauser (KS) 
function (ref. 10). 


Cij .= KS(gij) = i In [I exp (pg’j)] 

that has the property of approximating the maxi- 
mum constraint so that 


MAX ( g w ) < KS < MAX ( g w ) + ~ ln(W) 

with the factor p controlled by the user. Thus, 
the KS function serves as a convenient single 
measure of the degree of constraint violation (or 
satisfaction). 


In addition to the above functions the 
following functional relationship exist: 


XU 


Y ij 


Vector of the design variables, 
X’j, selected among the cross- 

sectional .. dimensions; t = 1 

through T 1 j . 

Vector containing the entries, 
Y’ J , of the stiffness matrix 

K 1 j and the mass M’j to be 
held constant; z = 1 through 
Z’J. 


g 1 j = f(X i j , Y 1 j , Q’j) 


The optimization problem definition is 


min C i j(X’j, Y’j, Qij) 
X^ 

h 1 j = 0 

L’j < X’j < U’J 


(13) 

(14) 

a) 

b) 

c) 
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In this problem, Y 1 J is being held constat by 
virtue of equation 14b, Consequently, Q’J re- 
main constant also and together with Y 1 J form 
a set of parameters of the optimization problem. 

Solution of this optimization problem (by 
any technique available) yields .a constrained 
optimum described by a vector n’J composed of 
th$ minimum value of the cumulative constraint, - 
C’J, and the optimal vector of the design 
variables, X 1 J 

n ijT = [cU|xU]T n .v 


includes^the constraints g”* representing the 
limits that may be imposed in the SSmkl own 
behavior, e.g., the interior displacements. It 
also includes the quantities t’J to account for 
the changes in C’J caused by the variability of 
the X"*. Ordinarily, every variation in X mk 
would require a reoptimization of the affected 
daughters SSijk. in order to find the new values 
of (TJ and X’J. However, following the con- 
cept introduced in reference 4 "‘and applied in 
reference 9 these new values will be approximated 
by a linear extrapolation using the optimum 
sensitivity derivatives. Thus, the potentially 
costly reoptimizations are bypassed. 


This solution i$ . sensitive to the parameters of 
the problem, Q 1 J and Y 1 J . That sensitivity, 
which will be used to approximate the daughter- 
parent coupling, may be quantified, as in refer- 
ence 9 by means of the optimum sensitivity 
derivatives (ref. 5). Considering that Q’J 
depends, on Y’J (equation 3 and the definition 
of X’J), the total derivative of C'J with 
respect to Y’J is 


Formalization of the above optimization 
problem follows the pattern established for the 
daughter substructure SSijk. 

Definitions: 

X mk Vector of the design variables, 

X™, t = 1 through T mk , 

related to the vectors Y’J by 


dc’jj s ac’J r &ll 

dY T J 3Y’J 1 3Q’J 3Y’ J 

z z r r z 


(16) 


In equation 16 the partials of C*,J with respect 
to Y^J and with respect to Q’^ are obtained 

from the algorithm described in reference 5, and 
the partial Q’ r J with respect to Y Z ’J can be 

calculated by conventional structural sensitivity 


analysis. Similarly, 



dX’J . ax’ J . v ax’j 
dY’J aY’I p3Q’ Jp 



3Y’J Z 

(17) 


Optimization Of A Parent Substructure 


As shown in figure 2, the parent substruc- 
ture SSmkl, m * 1-1, receives from its daughters, 
SSijk, the minimized values of their cumulative 
constraints, C’J, optimal values of their 
design variables, X’J, and the optimum sensi- 
tivity derivatives of thqs.e quantities with 
respect to parameters, Q’J and Y’J, accord- 
ing to equations 15 through 16. The substruc- 
ture SSmkl Itself is a daughter of SSnlp, n * m - 
1. that acts on SSmkl with the boundary forces, 
Q 1 *, and governs its boundary stiffness matrix 
by means of the parameters Y™. 


X mkT = [ ... j Y’J • ... ] T 


(18) 


Y”* Vector containing the entries, 

Y^ , z = 1 through Z"* , of 

the boundary stiffness matrix, 
K bmk , and the mass, M mk . 
These quantities are to be held 
constant in the ensuing optimi- 
zation. and in order to formulate 
an appropriate equality con- 
straint one needs to recognize 
the functional relationships 
discussed below. 


The matrix K bnik is a function of the 
stiffness matrix K™ 


K bmk , f ( ) 


(19) 


This function is no longer a simple identity as 
it was in equation 6, but represents elimination 
of the interior degrees of freedom by means of a 
solution of a set of simultaneous linear 

equations with many right-hand sides (e.g., 
refs. 6, 7, 8). 

In equation 19, the matrix K mk in turn is 
a function of the entries of the matrices K’J 
through a stiffness summation algorithm 


The design freedom in SSmkl consists in the 
freedom to manipulate the stiffness and mass 
distributions among its daughters by means of the 
design variables X^ collected in a vector 

X mk . That vector includes as partitions all 
the vectors Y’J whose entries were parameters 
in optimizations of the daughters SSijk. 


K mk = S(K’ J ) 

Hence, by equation 6 
K mk = f(X’J) 


( 20 ) 


( 21 ) 


Following equations 5 and 7, the mass, M mk , is 


The optimization problem in SSmkl is basi- 
cally the same as the one formulated previously 
for each daughter SSijk. It calls for finding a 
vector X"* that minimizes a cumulative con- 
straint, C mk , for SSmkl. That constraint 


M mk = f(Xij) = l M’j 

j 


( 22 ) 
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Q"* Vector of the boundary forces, 

Qmk. r = i through R"*. 

g"* Vector of the inequality con- 

straint functions, g™; w = 1 

through V/" k . 

h mk Vector of the equality con- 

straint functions, h™; s = 1 

through S"*. These constraints 
link Y"* and X"* so that, 
considering equations 18 through 
22 . 


h mk , 

ymk 

- f(K b ) * 0 

(23) 

hence 




h mk = 

Y mk 

- f(X mk ) = 0 

(24) 

L mk. 

U"* 

Vectors of 

lower and upper 


limits on the design variables 
Xj*. In this application, the 

L mk limits are needed tQ.keep 
the diagonal entries of lOJ and 
the masses M'J as nonzero, 
positive values, and also to 
prevent the off-diagonal entries 
of K’J from assuming physically 
unrealizable values. 


gink = f(x mk ( ym^ qmk) 


( 27 ) 


c;j • f(X«*. Y mk , Q"*) 


(28) 


Based on the above definitions, the optimization 
problem for the substructure SSmkl is 


mi n C mk ( X mk , Y"* , Q mk ) 
X mk 

h mk = 0 

L mk ^ x mk £ umk 

Lij £ )^lj £ uij 

where 

X ij = X 1j + A x mk 

e dX™ 

The increment aX’J is 
aX"* = X"* - Xo"* 


(29) 

a) 

b) 

c) 

d) 


(30) 


(31) 


duq to equation 30 in which X?* stands for 
Yp because of equation 18. 


Another Important role of the upper and 
lower limits is to represent the move limits 
needed to prevent excessive errors In the linear 
extrapolation of C'J and X'J. 

C e Vector of the daughter cumulat- 

ive constraints, t’J, estimated 
by a linear extrapolation and 
Included In the cumulative con- 
straint C . The extrapolation 
Is accomplished by a linear 
portion of the Taylor series 
using equation 16 taking Into 
account that the Xj* replaces 

the Y^ by virtue of the defi- 
nition in equation 18, so that: 


5ij . C 1 J + l AX"* 

e £ dX™< 


The constraints expressed by equation 30 are 
Introduced to reduce the probability that the 
optimization of SSmkl might Induce overstepping 
of the side constraints In the daughters SSijk. 
However, evaluating these constraints for the 
quantities approximated by equation 30 does not 
guarantee that such overstepping will not occur. 
Other means are needed to prevent that as it will 
be explained in the discussion of the entire 
iterative procedure. It follows that the con- 
straints of equation 29d are not essential and 
may be omitted. Indeed, there Is a strong motiv- 
ation to omit them because to do so would allow 
limiting the optimum sensitivity analysis to cal- 
culation of the objective function derivatives. 
This according to reference 5, requires an Input 
of only the behavior gradients while calculations 
of the optimal design variable derivatives 
requires Input of the gradients and the second 
derivatives of behavior - a very substantial 
different in the computational cost. 


C"* A cumulative constraint for 

SSmkl representing all functions 
in the vectors g"* and C e 
through the KS function intro- 
duced In equation 11, so that 


C mk = 1/p In [J exp( P g w ) 
w 


+ I exp ( pC 1 J ) ] 
ij e 


(26) 


Finally, owing to equation 24 there are 
functional relationships analogous to equation 13 
emended to include the cumulative constraints 

c ; j 


In view of the above, a distinction will be 
made between variant 1 and the algorithm that 
includes equation 29d and variant 2 in which It 
Is omitted. For completion of the presentation, 
the algorithm description will continue for 
variant 1. 

The SSmkl optimization produces a constrain- 
ed optimum described by a vector n"* composed 
of the minimum value of the cumulative con- 
straint, C"*, and the optimal vector of the 
design variables, X"*. 
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„mkT „ [5mk;jin*]T (32) 

It is followed by the analysis of sensitivity 
with respect to parameters Y mk and () m k con- 
sidering that the forces Q mk depend on Y mk 
through analysis of SSm!<l. The optimum sensi- 
tivity derivatives are 


5. There is no need for the equality con- 
straints to enforce constancy of the 
mass and the boundary stiffness entries. 

The definition needed to formulate the top 
level optimization problem are identical to those 
given for SSmkl, omitting those that do not apply 
because of the differences 1 through 5 above. 
The remaining definitions are:.. 


riC lllk 

3C" lk 

r 

3 C mk 

3() n,k 



XU 

Ve 

:ctor of 

the 

dY lTlk 


5 o*' 

almr 


( 33 ) 


*, 

11 t 

t ’ L 

= 

2 

Z 


r 

Z 



related to 

the 

dX' nk 

v'/inl: 

C A 


ax"* 

3Q mk 



X 11 T 

• l ... 

: y 2 j : .. 

• i 

. ] T 

Jy'mk 

z 

gyinF 

+ 1 
r 

aof 

a yml< 


( 34 ) 

The mass. 

MU i s 



which by 

vi rtue 

of 

recurs 

ivity of 

the 

decompo- 

MU ; 

= f(xii) 

= l M Z J 


sition are identica 

1 to 

equations 

16 

and 17 , 



j 



1 through T 


,2k' 


except the indexes ij being replaced by mk. 

The data defined by equations 32, 33 and 34 
are carried from SSrn!:l to its parent SSnlp at the 
levs 1 n - rp - 1 . 




( 35 ) 


( 36 ) 

Vector of the external loads, 
Qll, reduced to the boundary 


Opt i miza ti on 0 f T he Next Parent Structu re 

Moving on to the substructure SSnlp, every- 
thing that was stated in the preceding subsection 
on optimization of SSmkl applies to SSnlp 
literally, provided that: the indexes n, 1, and 
p are replaced by another triplet, say, a, p, c, 
that identifies the parent of SSnlp at the level 
a = n-1; and the indexes m, k, and 1 are replaced 
by n, 1, and p, For generality of variant 1, one 
needs also to extend equation 30 to encompass 
fully each line of succession emanating downward 
from SSnlp. Beyond these changes, no new concep- 
tual elements are introduced, and no additional 
definitions or discussion are needed at the junc- 
tions between the levels until one arrives at the 
top level. Hence, any number of intermediate 
levels of substructuring can be inserted, if 
physically justified, into a line of succession 
connecting the assembled structure on the top to 
any most detailed substructures below. This 
property characteri zes the algorithm as recur- 
sive. 


degrees of freedom according to 
the conventional formulation of 
s^structuring; r - 1 through 

gll Vector of the inequality con- 

straint functions, g 11 , im~ 

w 

posed on the displacements, 
buckling critical loads, or 
internal forces in SS 1 10; w = 1 
through wH. 

1.11, gll Definition for L mk , U"* that 
followed equation 24 applies with 
the indexes ij replaced by 2 j , 
and the indexes mk replaced by 
11 . 

C e Definition given in conjunction 

with equation 25 applies with the 
index replacements as above. 

Finally, based on the above definitions, the 
optimization problem at the top level is 


Op timization At The Highest Level 

* 

The assembled structure is designated 
SS110, Its optimization problem is similar to 
the one described for SSmkl with the following 
differences: 


min Mll(xH) 
XU 

gll < 0 


(37) 

a) 

b) 


1. No parameters are defined solely for the 
decomposition purposes. 

2. The objective function is the mass of 
the assembled structure. 

3. There is no need for a cumulative con- 
straint although there still is an 
option to use it to reduce the number of 
constraints that need to be processed. 

4. The boundary forces are the external 
loads on the assembled structure. 


Cfj< 0 c) 

L 11 £ yll £ yll d) 

l z i < 'x2j < U 2 J e) 

where equation 37e is analogous to equation 29d 
with the limits L‘J, U Z J reflecting the 
limits passed upwards through extrapolations of 
the type expressed by equation 30 extended 
recursively to encompass all the levels below as 
mentioned in the subsection on SSnlp. 
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Unlike all the daughters SSijk, the optimi- 
zation of SSI 10 does not have to be analysed for 
the optimum sensitivity. 

Iterative Procedure 

When the SSI 10 optimization is completed, 
the entire structure has acquired a new distri- 
bution of stiffness and mass within the move 
limits. Hence, the analysis must be repeated and 
followed by a new round of substructure optimi- 
zations in an iterative manner until conver- 
gence. The iterative procedure is composed of 
the following steps: 

1. Initialize all cross-sectional dimen- 
sions. 

2. Perform a substructuring analysis, 
including for each substructure at each 
level the transformation of the stiff- 
ness matrix into the boundary stiffness 
matrix, and the transformation of the 
forces applied to the interior degrees 
of freedom to the forces coinciding with 
the boundary degrees of freedom. 

Calculations of the behavior derivatives 
needed for the ensuing optimizations and 
for the optimum sensitivity analyses are 
implied in the substructuring analysis. 

3. Perform the operations of optimization 
and optimum sensitivity analysis as 
defined by equations 10 through 34. 

4. Optimize the assembled structure as 
defined by equation 35 through 37. 

5. Repeat step 2 .and terminate when: all 

constraints gf*J are satisfied at all 
levels, and M 11 has entered a phase of 
diminishing returns. Otherwise, 

continue. 

Salient Features Of The Algorithm 

Taking a perspective view on the multilevel 
algorithm with the algorithm for optimization 
without decomposition as a reference, the follow- 
ing salient features stand out: 

1. Optimization by decomposition replaces 
optimization of a single large problem 
with optimizations of a multitude of 
smaller problems which are isolated from 
each other. Obviously, all the sub- 
problems at a given level can be 
analyzed and optimized simultaneously. 

2. The design variables above the most 
detailed substructures are generalized 
design variables that enable the design- 
er to control the structure behavior by 
controlling its stiffness and mass 
distributions. 

3. The structure mass is controlled at the 
assembled structure level. The optimi- 
zations at the lower levels are concern- 
ed with improving the constraint satis- 
faction. 


4. The coupling among the subproblems is 
approximated by estimating the changes 
in the substructure behavior caused by 
the changes in the higher level sub- 
structures by means of linear extrapol- 
ations. In this respect, the approach 
resembles the piecewise linearization 
technique known to be effective in many 
optimization application (s.g., ref. 3). 

5. The multilevel and corresponding single 
level optimizations can be regarded as 
equivalent in the sense discussed in 
(ref. 9, App. B). Although no rigorous, 
optimality-condition based proof has as 
yet been produced, the equivalence 
assertion is supported by all the appli- 
cation experience to date. The equival- 
ence means that a multilevel optimi- 
zation and a single level optimization 
of the same problem will arrive at the 
same solution, if the problem is 
convex. In a non-convex problem, 
different solutions are likely because 
the two algorithms are different and can 
follow different search paths in the 
design space. 

6. From the view point of the overall opti- 
mization procedure, the operations of 
substructure analysis, optimization, and 
optimum sensitivity analysis are "black 
boxes" whose content can be freely 
replaced as long as their input/output 
remain as defined in the foregoing. 
Dissimilar algorithms may be used for 
the same operation applied to different 
substructures. In particular, although 
one would expect a finite element method 
to be used throughout for substructuring 
analysis, it is not a requirement as 
illustrated by the numerical example in 
the next section. 

7. The design freedom (equation 8 and 9) 
must exist between each parent and 
daughter substructures. Normally, that 
freedom is assured, if none of the 
daughter substructures is of the 
ultimate detail type. This is because 
the lack of such freedom would be an 
indication that the daughters do not 
contribute enough stiffness entries to 
the parent to support all of its degrees 
of freedom - a sign of incorrect sub- 
structuring. Given a design freedom, 
there still is the issue of the 
completeness of the control a designer 
wants to exert over the stiffness and 
mass distribution in a substructure. 
That control may be complete, if all 
entries of the substructure boundary 
stiffness matrix and the mass of a 
daughter substructure are design 
variables in the parent substructure. 
The control may be incomplete if one 
decides to manipulate as design 
variables only a subset of these 
quantities leaving the remainder free to 
float. The degree of control that one 
needs is probably problem-dependent and 
involves engineering judgment. 
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Intuitively, it appears likely that in 
many applications an incomplete control, 
e.g., control over only the diagonal 
entries of the boundary stiffness matrix 
and the mass, should be adequate. In 
some cases, a complete control may not 
be possible physically, e.g., one can 
not control all stiffnesses of a gener- 
ally anisotropic composite material 
plate, if there are not enough plies in 
its layup. The completeness of the 
stiffness and mass control remains to be 
a research issue to be investigated 
further. 

Numerical Example 

The algorithm presented in the preceding 
sections was tested on the portal framework 
structure that also served as a test case in 
reference 9. The structure is illustrated in 
Figure 3. In the original construction used in 
reference 9 the three beams of the framework had 
a thin-walled, I-shaped cross-section, therefore, 
the structure decomposed into two levels: 
assembled framework and individual beams. Figure 
3 shows how the construction was modified by 
replacing the I-cross-section beams with box 
beams. Each box beam is built-up of four thin, 
stringer reinforced walls in order to provide a 
third decomposition level. Although the struc- 
ture appears simple, the experience that has 
accumulated since it was introduced in reference 
9 shows (references 11 and 12) that it is a chal- 
lenging optimization test replete with local 
minima. 

Testing of the algorithm was carried out by 
first optimizing the structure without decompo- 
sition in a conventional- manner to establish a 
reference and then optimizing it as a three-level 
system. Each multilevel optimization was started 
from different points in the design space and the 
results were evaluated against the reference 
results. 

Description Of The Test Problem 

The construction, loading, design variables, 
constraints and the objective function of the 
test structure are: 

1. Construction: As shown in Figure 3. 

The box beams are symmetric with respect 
to the plane of the figure. The 
material is an Al-alloy with properties 
given in Appendix. 

2. Loading: Figure 3 shows the loads 

applied to the structure. The 
concentrated force and the moment 
constitute two, independent loading 
cases. 

3. The design variables are the cross- 
sectional dimensions labeled 
through X B g in Figure 3, DETAIL 8 and 
X(J through X^ in Section AA. The 
thicknesses 1, 3, and 4 In a box beam 
are equivalent wall thicknesses and 
incorporate the "smeared" stringers so 
that the thickness of the sheet metal 


itself is not an independent variable 
but results from the stringer dimensions 
and the equivalent wall thickness. The 
total number of the independent design 
variables is, considering symmetry of 
the beam cross-section: 6*3*3+2* 

3 = 60. The constraints are evaluated 
using the analysis tools described in 
the next subsection. 

4. Constraints: Horizontal translation and 

rotation at the loaded corner are limit- 
ed to .25 in and 0.005 rd, respectively. 
The beams must not buckle in a Euler 
column mode. The stiffened plates in 
the box beams must not develop stresses 
higher than the allowables stated for 
the material, and they must not fail by 
local buckling. There are also minimum 
gage constraints on the thicknesses and 
side constraints on the other cross- 
sectional dimensions. The cross- 
sectional geometrical proportions are 
also restricted by such obvious 
considerations as the need to keep each 
stringer from protruding too far into 
the interior of the box beam. The total 
number of constraints was 265. 

5. The structural material volume is the 
objective to be minimized. 

Analysis and Design Space Search Tools 

Calculation of displacements of the frame- 
work joints and the end-forces acting on its 
three beams was carried out using a small finite 
element program based on the displacement 
method. Each beam was modeled as a single beam 
element characterized by its cross-sectional area 
and bending moment of inertia in the framework's 
plane. In this phase of analysis the problem was 
treated as two-dimensional and each nodal point 
had only three elastic degrees of freedom. The 
support points were assumed clamped. The program 
was capable of calculating analytically the 
static behavior first and second order sensi- 
tivity derivatives with respect to the cross- 
sectional properties of area and moment of the 
inertia. 

Analysis of an individual beam under action 
of the end-forces was carried out by strength of 
materials relationships to obtain distributed 
normal and shear edge forces acting on the beam 
walls and to compute the corresponding average 
stresses. The beam column buckling was analyzed 
by a “designer handbook" type of closed form 
formulas provided in reference 13, based on the 
classical Eulerlan approach. 

Analysis of an Individual stiffened plate 
extracted from the box beam .included calculation 
of stress for the given edge-forces, equivalent 1 
Von Mises-Huber stress, and evaluation of several 
local buckling modes. These modes accounted for 
buckling of the sheet metal between the stringer 
and the plate edge, column buckling of the 
stringer with a cooperating strip of the sheet 
metal, buckling of the stringer web and flange, 
and flexural -twist buckling of the stringer 
elastically restrained by the sheet metal. A 
limited post-buckling analysis in the elastic 
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range was also included. The buckling analysis 
was carried out at the level characterized by 
reference 14 and it w’s coded in a program 
described in detail in ference 15. Oetailed 
information on the side constraints on the design 
variables and on the geometrical proportions of 
the cross-sections is provided in the Appendix. 

A general purpose program (ref. 16) based on 
the technique of usable-feasible directions was 
used for the design space search in the reference 
optimization and at every level of the three- 
level optimization. 

Three-Level Optimization 

To establish the reference results, the • 
framework was first optimized without decompo- 
sition. Then, the multilevel optimization algor- 
ithm, variant 2 was applied to the structure 
decomposed as shown in Figures 3 and 4 showing 
the stiffened panels as daughters clustered in 
triplets (the fourth wall is symmetric) under a 
parent box beam. The beams, in turn, are daugh- 
ters of the assembled structure. 

The definition of the objective function, 
design variables, and constraints for each level 
in the decomposition is given in Table 1. As 
shown in the table, the top level optimization 
manipulates the beam extensional and bending 
stiffnesses through the cross-sectional area and 
bending moment of inertia. By coincidence, the 
area controls also the beam volume which contrib- 
utes directly to the objective function. 

At the middle level, the stiffnesses 
expressed by the area and moment of inertia 
become fixed parameters and the variables are the 
wall membrane stiffnesses controlled by the 
geometrical dimension variables. These vari- 
ables, and consequently the membrane stiffnesses 
become fixed parameters at the bottom level at 
which the ultimate detail dimensions are engaged 
as variables. The equality constraints arise 
between the parameters and variables. Owing to 
relative simplicity of the expressions involved, 
(see Appendix) these constraints are solved 
explicitly. 

Examination of Table 1 in conjunction with 
the previous description of the analysis tools 
illustrates the point that dissimilar analysis 
may be used as needed at different places in a 
decomposition scheme. 

The sensitivity analysis of behavior has 
been carried out by a single step forward finite 
difference technique. The optimum sensitivity 

analysis was based on the algorithm given in 
reference 5. 

Results And Remarks On The Method Performance 

Figures 5, 6, and 7 show a sample of results 
obtained when starting from with and without 
decomposition. To assure comparability of the 
results the starting points for both methods are 
the same. The normalized plots illustrate for 
each of the three different starting points, the 
objective function, a selected individual 


constraint, and a cumulative constraint contain- 
ing the above individual constraint as they 
varied over the iterations in the optimization 
without decomposition and cycles in the three- 
level optimization. An iteration is defined in 
the optimization without decomposition as a 
usable-feasible directions iteration. A cycle is 
defined as one execution of the series of steps 
listed in the iterative procedure definition in 
the previous section. 

The results verified that the multilevel 
algorithm was capable of finding a feasible 
design having an objective function close to and, 
in some cases, lower then the reference optimi- 
zation without decomposition. As in reference 9, 
differences up to 72.1% were observed among the 
detailed design variables obtained by the two 
methods. However, these differences were no 
larger than those observed by comparing the 
designs obtained without decomposition starting 
from different initial design points. Therefore, 
these differences can be attributed to the 
problem non-convexity. 

The volume of the data needed to describe 
the optimization results for the test structure 
is so voluminous that only a sample for one opti- 
mized case corresponding to Figure 6 is displayed 
in Table 2 to show dimensions of the optimized 
cross-sections. 

The objective function minimum in the 
three-level optimization falling below the value 
obtained by the optimization without the decompo- 
sition was an unexpected result that occurred in 
several tests two of which are illustrated in 
Figures 5 and 6. Examination of the detailed 
numerical data suggested that these particular 
results were caused by both the different search 
path taken and by the larger number of con- 
straints that the usable-feasible directions 
search algorithm had to process simultaneously in 
the optimization without decomposition. The 
algorithm implemented as described in reference 
16 does not have a rigorous, Kuhn-Tucker based 
termination criterion. Instead, it terminates by 
a "practical 1 ' criterion of the diminishing 
returns in the value of the objective function. 
The criterion numerical tolerance was set the 
same for both methods. Proceeding from one con- 
straint boundary to the next, as the usable- 
feasible directions algorithm does, the objective 
function reductions measured between the consecu- 
tive iterations tends to fall below that 
criterion prematurely, if there are many closely 
spaced constraint boundaries. In contrast, the 
multilevel scheme incorporates a piecewise 
linearization that inherently tends to produce 
significant increments of the objective function 
from one cycle to the next and tend to proceed 
further toward the theoretical constrained 
optimum than the reference method. * 

The graphs in Figures 5, 6, and 7 have a 
jagged appearance for both methods which is a 
characteristics of the usable-feasible directions 
search algorithm, amplified in the multilevel 
optimization by the extrapolation errors. How- 
ever, these errors have never become excessive 
and the daughter substructure reactions to the 
changes in the parent design were effectively 
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predicted by the optimum sensitivity deriva- 
tives. In at least one case these predictions 
enabled the optimization at the middle level to 
remove the constraint violations at a bottom 
level substructure without any change to that 
substructure sizing. 

Regarding the computational efficiency, it 
was not the purpose to demonstrate improvements 
of that efficiency and neither the reference nor 
the multilevel optimization procedures were honed 
for best computational performance in their 
implementations. The test case was compu- 
tationally too small anyway to permit drawing 
conclusions as to the computational efficiency of 
any method. The only observation in this regard 
is that the amount of computational labor in one 
cycle tends to be less than in one iteration. 
Therefore, the total numbers of cycles and iter- 
ations should not be compared to evaluate the 
multilevel method efficiency. The main advantage 
of the multilevel method stems from its compati- 
bility with the distributed computing technology. 

Concl usions 

An t algorithm has been described for perform- 
ing structural optimization by decomposing an 
optimization problem for an entire structure into 
a set of smaller subproblems. Each subproblem 
corresponds to a substructure in a general, sub- 
structuring analysis based on many levels of 
nested substructures. The optimization sub- 
problems remain coupled by means of the optimum 
sensitivity analysis that generates derivatives 
of optimum solution with respect to constant 
parameters. These derivatives quantify the rate 
of the behavior change in a substructure opti- 
mized for constant inputs received from its . 

governing, higher order substructure relative to 
the rate of change of these inputs. This quanti- 
fication informs the optimizer modifying the 
higher order substructure about the effect its 
action will have on the subordinated substruc- 
tures. The algorithm is ihtrinsically compatible 
with distributed computing because the sub- 
problems are isolated and can be processed in 
parallel . 

Testing on a framework structure made up of 
box beams that decomposes into a three-level 
pyramid of 13 subproblems showed that the optimi- 
zation algorithm performed as expected or 
better. Such difficulties as have been encoun- 
tered were not of the computational nature but 
were caused by complexity of data handling by the 
conventional FORTRAN means of named files and 
COMMON blocks. In fact, from the programmer's 
viewpoint this was a data-doml nated problem and a 
conclusion was drawn that to proceed with a 
larger scale applications, with more levels and, 
particularly, with more engineering disciplines 
involved in the analysis would require systematic 
means of data handling such as those provided by 
modern data base management systems (e.g., ref. 

17). 

Comparison against a reference optimization 
without decomposition of the same structure 
demonstrated the validity and effectiveness of 
the multilevel optimization algorithm. 


Satisfactory testing of the three-level 
optimization is seen as verifying an entirely 
general ,• multilevel algorithm because the formu- 
lation is recursive, and extension beyond three 
levels introduces no conceptually and qualitat- 
ively new elements into the decomposition 
scheme. The reported implementation is, there- 
fore, regarded as a stage in the development of a 
multilevel, multidisciplinary optimization scheme 
applicable to engineering systems such as air- 
craft, due to the generality of the basic con- 
cepts underlying the algorithm. 
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APPENDIX: Details- of the Formulation 
Additional or Redefined Nomenclature 
A z Z-stiffener cross-sectional area 

a Z-stiffener leg cross-sectional area 

d X2 or Xg 

Iy Moment of inertia about the axis 

parallel to the webs 

L Move limit percentage 

i Beam length 

M Bending moment 

N Resultant axial force 

P Axial force 

Q First moment of the area 

R Thickness-to-length ratio of the 

z-stiffeners legs 

S Section modulus 

s Slenderness (equation A10) 


t Smeared thickness 

V • Shearing force 

y Centroidal distance to the extreme fiber 

a Normal stress 

t Shearing stress 

Miscellaneous Information 

The material was 6062-T62 aluminum defined 
by Table 48, reference (13). The box beams were 
assumed to be singly symmetric, allowing only one 
web to be modelled. The cross-sectional area of 
the z-stiffener was assumed to be a constant 10% 
(r = 0.9) of the smeared area of the flanges and 
5% (r = 0.95) of the webs. The panels were 
assumed to be braced for local buckling every 24 
Inches. 

Details of the Optimization Without Decomposition 

Six design variables were used to model each 
z-stiffener. Six more design variables were used 

M M 

for the beam widths, X 2 , and heights, X5. The 
other beam dimensions were calculated by 

A z = X*X B 2 + X®X| + X B (X® - X B - X B ) A1 

t = A z /[b(1.0 - r)] A2 

M M 

where b * X 2 for flange panels or X5 for web 
panels. Beam section properties were calculated 
from the dimensions. A finite element program 
was used to analyze the structure for displace- 
ments and end forces for the two independent load 
cases, 50,000 lbs. and 20 x 10® lb-1 n. 

Four displacement constraints were computed. 


for 

each load case: 




D J 

c j 

= 5 for each load 
■ 0.25 in. 

case, 

j - 1,2 

A3 

C J 

= e for each load 
a 0.005 rad. 

case. 

j = 1.2 

A4 

9i 

3 Oj/cj - 1.0 

i = 

1,4 

A5 


For both loading cases, the stresses were 
evaluated at six points. Bending stresses were 
calculated at the bottom and top of both ends and 
shear stresses were evaluated at the neutral axis 
at both ends. 



When evaluating the stress constraints, a factor 
of safety was appl i ed 

Dj * l.‘5o 

Cj * 26,000 psi for tension 
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Dj = 1.5 


15,000 psi 


Dj = X?, x®. X® for webs 
A9 Cj = (X 2 - 2X^)/4 


For compression, the allowable stress was reduced 
to prevent buckling. The compression allowable 
was computed by 

s 3 S1/I y A10 


Dj = X i , X®, X® for flange 

Cj - X5/4 


A21 


A22 


For s < 40 Cj = 26,000 psi 

For 40 < s < 1810 Cj ■ 1,000(27.4 - 0.22/I)psi 

For s < 1810 cj = 1810 (18040)/s psi AU 

The function for the compression allowable was 
patterned after Table 4b of the ASCE specifi- 
cation for 6062-T6 aluminum in reference 13. The 
36 stress constraints were 


Dj= xjj, X®. X® for webs 
Ci = ^/4 


gi = Dj/Cj - 1.0 i = 212, 265 


A23 

A24 


Upper and lower- bounds on the design vari- 
ables were 


gi = Oj/Cj - 1.0 i * 5, 40 


A12 


5.0 < X 2 , X5 < 100.0 for all cases 


A25 


The end moments were converted to couples 
and summed with the axial force by 


N{ - Pi + M 1 /(X I J + X^ + X?) 


A13 


The resultant forces on both ends of the beam for 
both load cases were applied as four separate 
loadings. The shears and resultant forces for 
both flanges were applied to the webs, generating 
eight loadings. Thirteen constraints for each 
load case were evaluated to select Wie thirteen 
maximum values. This approach resulted in 117 
constraints for all nine panels. 

Additional geometric constraints were 
applied to the z-stiffener in the following forms 


a i * X®X®, a 2 s X3X5, a 3 * X®(X® - X® - X®) 

A14 


Dj * a i 
Cj " 0.1A Z 


gi = Dj/Cj - 1 i = 158, 184 


A15 

A16 


Ri = X®/X®, R 2 « X®/X®, R, = xf/xf 


A17 


°j * R X 

cj = 0.5 


A18 


9i * Dj/Cj - 1.0 i » 185, 211 


A19 


The z-stiffener dimensions were also constrained 
to fit Inside the box beam by 

D. * X®, X®, X® for flanges 

J 


Cj - (X? - 2 x ”)/4 


A20 


0.0625 < X®, x|, X® < 2.5 

1.0 < X®, X® < 25.0 

1.125 < X® < 25.0 


for case 1 
A26 


0.0625 < X®, X®, X® < 
1.0 < X®, X® < 

1.125 < X® < 


1.0 


10.0 for case 2 and 3 


10.0 


A27 


The objective function, F, was defined as 
the framework material volume. 


F 


3 


= I A i 1 i 
i=l 


A28 


Details of the Multilevel Optimization 


Details of the bottom level formulation have 
been given in the preceding section. At the 
middle level if and U-^ were defined for X 2 
and ^5 by equation (A26) . For xj 1 , X^ and X^, 

BL^ and U 3 M were calculated by replacing X 3 

with U® and L® in equation (Al) to calculate 
Az, max and A z, min. respectively. Then 

b max( 1 - r ) A29 


Ui » 

bmin^" r ) 


A30 


where b m ^ = 5 and b max = 100 consistent with 


wnere o m1 n - a 
equation (A25). 


The more restrictive of U M and l M or the 
move limits, calculated by 
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A31 


X U,L = x M (l ± L) 
L = 0.075 


A3 2 


were chosen as the cycle's upper and lower 
bounds. 

MM M 

x", X 2 , and X 3 were used as design variables 
while X 4 and X 3 were calculated as dependent 
variables to satisfy equation (15) by 


B, = 0.5XM ) 2 - (X^)2]/A 


B 2 * 0.5X^[X^ - Xi]/A 

B 3 = [A - X 2 (X^ + Xj)]/12 + X^X 2 (B 2 + 1/2 ) 2 
+ X^X^(B 2 - 1/2 ) 2 + B|(A - X^X 2 - X 2 X 3 ) 


A33 

A34 

A35 


B 4 = 2[X I Jx 2 (B 1 + X^/2) (B 2 + 1/2) 
+ X 2 X 2 (B! - Xj/2) (B 2 - 1/2) 

+ 8 38 2 ( A - X^X 2 - X 2 Xj)] 


A3o 


B 5 = X^X 2 (B j + Xj/2 ) 2 + XjX 3(8 j - Xj/2 ) 2 

+ B f ( A - X?X5 - X 2 X 3 ) + X 2 [(X ^) 3 + (X 3 ) 3 ]/ 1 2 - I 

A3 7 


X 5 = 0.5(-B 4 + - 4B 385)7 


/ B, 


X 4 * 0.5( A - X^X« - X M 2 X M 3 )/X M 5 


A38 


A39 


At the top level, the more limits were defined as 
in equations A31 and A32. 


Table 1. Quantities Defined for the Multilevel Test Case Optimization 



TOP LEVEL 

OBJECTIVE: 

The framework material volume 

DESIGN VARIABLES: A and I of the -beams 

CONSTRAINTS 

Displacements of the loaded corner and for the beams 


MIDDLE LEVEL 

OBJECTIVE: 

Cumulative constraint C representing the column buckling 

and C for the walls, 
e 

DESIGN VARIABLES: Wall membrane stiffness contributing to the beam 

axial and bending stiffnesses controlled through 
the dimensions shown in fig. 3, section A-A. 

CONSTRAINTS 

Equality-beam cross-sectional area and moment of inertia 


BOTTOM LEVEL 

OBJECTIVE: 

Cumulative constraint C representing a set of stress and 
local buckling constraints of the wall. 

DESIGN VARIABLES : Cross-sectional dimensions shown in fig. 3 DETAIL B 

CONSTRAINTS 

Inequality - minimum gages, geometrical proportions, and 
geometrical realizability. 


Equality - membrane stiffnesses for tension-compression 
and bending in of the wall in its own plane. 
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Table 2. REPRESENTATIVE RESULTS 



BEAM 1 

BEAM 2 

BEAM 3 



SINGLE 

MULTI- 



MULTI- 



MULTI- 


INITIAL 

LEVEL 

LEVEL 

INITIAL 

HI* Mil 

LEVEL 

INITIAL 


LEVEL 


VALUES 

FINAL 

FINAL 

VALUES 

FINAL 

FINAL 

VALUES 

FINAL 

FINAL 



VALUES 

VALUES 


VALUES 

VALUES 


VALUES 

VALUES 

F(in3) 

64347.5 

57127.0 

46060.5 

64347.5 

57127.0 

46060.5 

64347.5 

57127.0 

46060.5 

A(in i! ) 

42.7 

42.1 

36.3 

190.0 

161.8 

127.7 

57.7 

55.2 

46.1 

Kin 4 ) 

3903.5 

3894.0 

3365.0 

31088.1 

41370.0 

39550.0 

5528.1 

6169.0 

4621.0 

q 

y max 

0.092 

-0.097 

0.044 

0.092 

-0.097 

0.0438 

0.0917 

-0.097 

0.044 

C M 

- 


0.006 

- 

- 

-0.157 

- 

- 

-0.067 

Xj 1 ( 1 n) 

1.620 


1.572 

0.971 


0.816 

■ 0.791 


0.650 

(in) 

10.0 

9.997 

8.524 

35.0 

34.34 

31.32 

20.0 

19.96 

18.57 

X 3 ( i r>) 

0.75 r 


0.831 

1.3 


1.456 

0.88 


0.787 

X 4 (in) 

0.423 


0.354 

1.727 


0.677 

0.518 


0.420 

X^(in) 

22.5 

22.65 

22.35 

32.0 

38.71 

41.74 

22.5 

24.15 

23.08 

C 8 

- 

- 

0.024 

- 

- 

-0.182 

- 

- 

-0.041 

X?(in)' 

0.1 

0.1 

0.089 

0.8 

0.845 

0.851 

0.3 

0.297 

0.383 

X 8 (in) 

0.1 

0.1 

0.086 

0.25 

0.263 

0.325 

0.4 

0.396 

0.176 

X 3 (in) 

0.4 

0.398 

0.358 

0.8 

0.845 

0.851 

0.3 

0.297 

0.383 

X 8 (in) 

1.0 

1.0 

1.0 

2.0 

2.130 

2.176 

2.0 

1.976 

1.627 

X 5 ( i n ) 

1.0 

1.0 

1.004 

2.0 

2. 130 

2.176 

2.0 

1.976 

1.627 

xjdn) 

3.0 

2.988 

2.941 

7.0 

7.497 

2.176 

2.0 

1.976 

1.627 

c 8 

- 

- 

0.041 

- 

- 

-0.227 

- 

- 

-0.040 

X? (in) 

0.3 

0.299 

0.35 

0.5 

0.495 

0.285 

0.175 

0.174 

0.072 

xf (In) 

0.3 

0.296 

0.160 

0.4 

0.371 

0.305 

0.175 

0.174 

0.092 

X 3 (in) 

0.3 

0.299 

0.321 

0.5 

0.495 

0.165 

0.5 

0.491 

0.491 

X 8 (in) 

1.0 

1.0 

1.272 

1.0 

1.0 

1.0 

2.0 

1.987 

1.804 

X|(1n) 

1.0 

1.0 

1.240 

1.0 

1.0 

1.721 

2.0 

1.963 

1.667 

x|(in) 

4.0 

3.934 

3.420 

7.0 

6.404 

7.755 

2.0 

1.987 

2.523 

C 8 

- 

H 

- 0.110 

- 

- 

-0.159 

- 

- 

-0.144 

Xj(in) 

0.062 


0.086 

0.3 

0.285 

0.176 

0.062 

0.064 

0.062 

x|( 1 n) 

0.1 


0.063 

0.625 

0.513 

0.172 

0.062 

0.064 

0.069 

Xj(1n) 

0.175 

0.173 

0.155 

0.2 

0.187 

0.172 

0.35 

0.338 

0.229 

x 5 ( 1 n) 

2 . 0 ' 

1.986 

1.545 

2.0 

1.856 

1.940 

1.0 

1.0 

1.0 

xf (i n) 

1.0 

1.0 

1.004 

3.0 

2.784 

2.962 

1.0 

1.0 

1.016 

X 8 (1n) 

2.0 

1.977 


3.0 

2.343 

2.851 

3.0 

2.938 

3.131 

Iterations 

7 



7 

22-23 


7 

22-23 
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Level I 



Fig. 1 Multilevel substructuring. 
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Fig. 5 Representative results. 
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Fig. 3 A portal framework. 
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O Single level results - iterations 
□ Multilevel results - cycles 
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O Single level individueal constraint <44 

□ Multilevel cummulative constraint containig constraint #44 

® Single level individual constraint #83 

69 Multilevel cummulative constraint containing constraint #83 
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Fig. 6 Representative results. 
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O Single level results - iterations 
^Multilevel results - cycles 


O Single level individual constraint #83 
□ Multilevel cummulative constraint containing constraint #83 
® Single level displacement constraint 
BMultilevel displacement constraint 
0 1 
0 

-01 


Constraint value 
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Fig. 7 Representative results. 
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